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Abstract 



where ri , r% £ C(x, t) and C is an algebraically closed field of characteristic zero. 

We work in the setting of parameterized Picard-Vessiot theory, which attaches a linear differential 
algebraic group to such an equation, that is, a group of invertible matrices whose entries satisfy a 

system of polynomial differential equations, with respect to the derivation in the parameter-space. We 
will compute the ^-differential-polynomial equations that define the corresponding parameterized 
Picard-Vessiot group as a differential algebraic subgroup of GL2. 
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0. Introduction 

0.1 Background 

In differential Galois theory, one studies a differential equation with coefficients in a given differential field K 
(§1.1), by investigating the differential field extension of K generated by the solutions for the equation, together 
with their derivatives. This is in analogy with the classical Galois theory of polynomial equations, where one 
studies such equations by studying the splitting field of the polynomial. In both cases, the algebraic relations 
amongst the solutions for the equation are reflected in the algebraic structure of the group of automorphisms 
of the field extension over the base field. The differential Galois theory of linear differential equations was 
developed by Kolchin [18,20], putting earlier work of Picard and Vessiot on a firm modern footing. The differ- 
ential Galois group corresponding to a linear differential equation is a linear algebraic group over the subfield of 
constants of K (§1.1). 

In [3], a parameterized Picard- Vessiot theory for linear differential equations with parameters is developed in 
close analogy with Kolchin's differential Galois theory. This theory is a special case of an earlier generalization 
of Kolchin's theory, developed in [23]. In the parameterized theory, the base field K is a differential field, where 
there are now two kinds of derivation: principal derivations and parametric derivations, so that the subfield of 
constants, with respect to the principal derivations, is a differential field with respect to the parametric deriva- 
tions. The PPV field (Definition 1.8) of a system of linear differential equations, with respect to the principal 
derivations, is generated over the base field by the solutions for the system, together with their derivatives with 
respect to all (principal, as well as parametric) derivations. The PPV groups (Definition 1.9) which arise in 
this theory are linear differential algebraic groups (Definition 1.2), and their structure reflects the differential- 
algebraic relations amongst the solutions, with respect to the parametric derivations. Linear differential algebraic 
groups are the differential-algebraic analogues of linear algebraic groups — that is, they are subgroups of GL„ 
which are defined by the vanishing of systems of polynomial differential equations in the matrix entries. 

The parametrized Picard- Vessiot theory has been finding a diverse number of applications in other areas. We 
shall list a few examples, although we will not attempt to be exhaustive. In [12, §5], the theory is related to Gauss- 
Manin connections [17,24], leading to many potential applications in algebraic geometry. In [3, 1 1, 12, 25] the 
theory is applied to the study of isomonodromy, and in [26] it is applied to monodromy-evolving deformations, 
and to the study of some equations in mathematical physics; see [25, §1] for more references on this. In [13], 
the authors apply a generalization of this theory to the study of difference equations, which has numerous 
applications. See also [11], where the authors recast the theory in a new light and suggest several applications. 

0.2 This work 

We develop algorithms to compute the PPV group (Definition 1.9) associated to a second order homogeneous 
linear differential equation defined over C(x,t); that is, an equation of the form 

d 2 Y dY 

^ T + r 1 — +r 2 Y = 0, (3.1) 

where r\ , r 2 £ C(x,t), and the principal and parametric derivations are and 4, respectively. The key ingredient 
in our approach is the Reduction Lemma 2.1, which reduces the determination of the PPV group to solving 
some systems of linear equations. Our algorithms are Maple-ready, and we are currently working on writing a 
computer implementation. 

For technical reasons, 1 it has become a tradition in the theory to assume that the field of constants, with 
respect to the principal derivations, is differentially closed (§1.1). Therefore, we work over the field K := Kq(x), 
where Kq is a ^-differentially closed field extension of C(t). Working with differentially closed fields tends to be 
an obstacle in practical applications of the theory, since such fields are usually unnaturally large (see [35, §1]). 

'in [3], the authors impose this condition in order to ensure the existence of the PPV groups and extensions (see § 1 .4). 
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The groups produced by our algorithm are defined over C(t) (Theorem 4.2), and we only perform computations 
over finite algebraic extensions of C(t) — the differential closure plays a marginal role in our arguments. This is 
in agreement with the stronger and more general results of [1 1, Theorem 2.8], which imply in particular that the 
groups we wish to compute are actually defined over C(t), while it is shown in [35], by different methods, that 
they are defined over C(t). 

In [22], Kovacic developed an algorithm which computes all the "elementary" solutions to a second order 
order equation of the form (3.2), if they exist. From these data, together with the classification of the algebraic 
subgroups of SL2, Kovacic's algorithm decides which algebraic subgroup of SL2 is the (non-parameterized) 
PV group (Definiton 1.7) of (3.2). The differential algebraic subgroups of Sh2(U), where U is a universal 
differential field, 2 were classified in [32]. In [3, p. 137], the authors ask whether one can use this classification to 
develop a parameterized analogue of Kovacic's algorithm. Although we do not rely on the classification itself, 
the methods in [32] were a great source of inspiration in developing the algorithms presented here, and the 
outputs of our Algorithms 4 and 5 have the form of one of the groups in the classification in [32] (Theorem 4.1). 

An algorithm to compute the PPV group of (3.2), provided that it is a Zariski-dense subgroup of SL2, is 
given in [8]. We develop algorithms to compute the PPV group of (3.1) in the remaining cases. Some of the 
arguments we give in the other cases are similar to those in [8], but we work in different settings: in [8], the field 
of principal constants is assumed to be universal, and a theorem of Seidenberg is applied in order to obtain a 
concrete analytic interpretation for elements of this abstract universal differential field, as meromorphic functions 
on some complex polydisk. Aided by this concrete interpretation, the author obtains his results. We will reprove 
some of these results algebraically in the course of the proofs of Propositions 3.5 and 3.7, in order to be able to 
apply them in our setting. 

Lastly, we wish to point to a subtlety which is new to the parametrized setting. In Kovacic's algorithm, one 
performs a change of variables on (3.1) in order to obtain a new equation (3.2), whose PV group is a subgroup 
of SL2. Then, it is a simple matter to reconstruct the PV group of the original equation from these data, which 
is always an almost direct product 3 of the PV group of (3.2) and a change-of-variables group (§3.4). This is 
no longer true 4 for parameterized equations, where the richer differential-algebraic structure of the differential 
Galois groups allows them to be related to each other in more complicated ways. We have solved the problem of 
determining precisely how the PPV group of (3.1) is related to that of (3.2) in all cases, and we give an algorithm 
to compute this relationship in §3.4. 

0.3 Organization 

We will now describe the contents of the paper in some more detail. In § 1 , we will briefly recall some basic 
definitions and results from differential algebra, the theory of linear differential algebraic groups, Picard-Vessiot 
theory, and parameterized Picard-Vessiot theory. In this section, we will also define most of the notation and 
concepts which will be used later in the paper. In §2, we will describe three auxiliary algorithms which will form 
the building blocks of our algorithms to compute the PPV group (Definition 1 .9). Each algorithm in this section 
is first followed by a concrete example of its implementation, and then by a proof that it gives the right answer. 
In §3, we will describe algorithms to compute the PPV group of (3.1). We will first apply Algorithms 1 and 2 
of §2 to compute the PPV group of an associated second order (3.2), which is a differential algebraic subgroup 
of SL2. Then, we will indicate how to recover the differential algebraic subgroup of GL2 corresponding to the 
original (3.1), as an application of Algorithm 3 from §2. In §4, we state some consequences of the algorithms 
and comment on improvements which will be relevant for practical applications. 



2 The notion of universal differential field [20] is stronger than that of being differentially closed. 
3 This is a quotient of the direct product of two groups, with finite kernel. 

4 This answers a question of Ovchinnikov [28], who suggested to the author to treat the general case where the PPV group is not assumed 
to be unimodular. He has pointed out to us that the parameterized and non-parameterized situations are still analogous: the PPV group 
of (3.1) is a quotient of the direct product of the PPV group of (3.2) and the change-of-variables group, with zero-dimensional kernel. 
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1. Preliminaries 

In this section, we will briefly recall without proof some standard definitions and results. We will also set the 
notation which we shall use for the rest of the paper. The algebraic closure of a field K will be denoted by K. 
References will be given in each section. 

Every field considered in this work will be of characteristic zero, 
and C will always denote an algebraically closed field. 

1.1 Differential algebra 

We will briefly recall some definitions from differential algebra. We refer to [16, 20, 29] for what follows. 

A differential field is a pair (K, A), where K is a field, and A is a finite set of pairwise commuting derivations. 
We abbreviate this by saying that K is a A-field. If A = {8} is a singleton, we write 8 instead of {8}, e.g. we say 
K is a 8-field. We denote the rc-fold composition of 8 with itself by 8", with the convention that 8 := id# . 

We say that a field automorphism a : K — > K is a A-automorphism if o o 8 = 8 o a for every 8 G A. For any 
subset II C A, we denote by K u the subset of elements c £ K such that 8c = for all 8 G IT. One can show that 
K is a field, and we call it the subfield of Yl-constants of K. 

If L is a A-field, and K C L is a subfield such that 8(K) C K for each 8 G A, we say that L is a A-field extension 
of K, and that K is a A-subfield of L. If y\, . . . , y m G L are such that 

is a set of generators for L as an overfield of K, we write L = K(y\ , . . . ,y m )A- In this case, we also say that L is 
differentially generated by {yi, . . . ,y m } over K, and that this is a set of differential generators for L over /f. 

Now let K be a 8-field. The ring of differential operators over K is denoted by K[5\ and consists of elements 
of the form Y,'"=o r n$ n with r n G K. The multiplication in this ring is given by composition, and is determined by 
the rule 8o/ = /o8 + 8(f) for / G K. The inte ger m is the order of the differential operator (D : — Y^n=o r n8 n , 
and we denote it by orA('D). We say T) is monic if r m = 1. 

The ring of differential polynomials over K (in m differential indeterminates) is 

K{Y h ...,Y m } s :=K[Y^}, 

the free (commutative) /f-algebra on the infinite set of variables {1^- | 1 ^ i ^ m, j £ 2^o}- We give it a 
8-ring structure by setting 8^ = Yp . We often omit the zero superscript, and write Yj = f/°' > . The ring 
^{Fi , . . . , F m }s is a free object in the category of 5-K- algebras. In particular, any ordered m-tuple (yi , . . .y m ) of 
elements of K defines a unique 8-homomorphism K{Y\ F m }s — > K, determined by Y^ \-t 8 ; y,. We say that 
K is differentially closed if every consistent system of differential polynomial equations with coefficients in K 
has a solution in K. 

1.2 Linear Differential Algebraic Groups 

The theory of differential algebraic groups was initiated by Cassidy in [1], and we refer to that paper for more 
information and complete proofs of the statements given here (see also [21]). For the rest of this section, we 
assume that K is a differentially closed 3-field, with subfield of constants K d = C. 

Definition 1.1. Let V C K m . We say V is Kolchin-closed if there are differential polynomials p\...p r in 

K{Yi,...,Y m } d such that „ , , N 

V = {v G K m | pi(y) = 0, V 1 < i < r}. 

Definition 1.2 (Cassidy [1, II, §5]). A linear differential algebraic group is a Kolchin-closed subgroup of 
GL n (K) cK nxn . 
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Theorem 1.3 (Cassidy [1, Prop. 1 1]). LetB be a proper differential algebraic subgroup ofG a (K). Then, there 
exists a unique nonzero monic operator L G K[d] such that 



Theorem 1.4 (Cassidy [1, Prop. 31 and its Corollary]). Let A be a proper differential algebraic subgroup of 
G m (K). Then, either A = p n C C x , the group of n th roots of unity, for some n G N, or else there exists a unique 
nonzero monic operator L € K[d] such that 



Theorem 1.5 (Cassidy [1, Prop. 42]). Let H be a Zariski-dense differential algebraic subgroup ofSL2(K). 
Then, either H = SLq,(K), or else H is conjugate to SL-2(C). 

1.3 Picard-Vessiot theory 

We refer to [29, Ch. 1] for what follows. Let K be a 8-field such that K 5 is algebraically closed and let CD be a 
differential operator over K of order m. 

Definition 1.6. A Picard-Vessiot extension (or FY -extension) for the differential equation CDY = is a 8-field 
extension of K, which we denote by PV(CD/K), such that 

(i) There are m distinct, /f 8 -linearly independent elements yi,...,y m £ PV(^D/K) such that CDy\ = for all 



A PV extension for the operator CD exists and is unique up to (non-unique) differential isomorphism, thus 
justifying the notation PY(CD/K). 

Definition 1.7. The Picard-Vessiot group of CD over K (or PY-group) is the set of 8-field automorphisms of 
PY(CD) which leave K point- wise fixed, and is denoted by Gal§(CD/K). 

The solution space of CD is the /f 8 -vector subspace of PV(CD/K) consisting of all y G PV(CD/K) such that 
CDy = 0, and we denote it by Sol(2>). 

The vector space Sol(©) is m-dimensional, and it is stable under the action of Gal§(©/,K). This action 
defines an injective homomorphism Gal$(CD/K) GL(Sol(£>)), and one can show that its image is a linear 
algebraic group. There is a Galois correspondence between Zariski-closed subgroups of Ga\$(CD/K) and inter- 
mediate differential field extensions between K and PV((D/K) (cf. Theorem 1.10). 

1.4 Parameterized Picard-Vessiot theory 

We are now ready to state the main results of the parameterized Picard-Vessiot theory. We will follow the 
presentation in [3]. Let K be a A := {S,3}-field, such that '■= is a differentially closed 3-field, and let 
CD G K[S\ be a differential operator over K, of order m. 

Definition 1.8. A parameterized Picard-Vessiot extension (or PFV 'extension) for the differential equation 
CDY = is a A-field extension of K, which we denote PPV(f /K), such that 

(i) There are m distinct, ^o-linearly independent elements yi,... ,y m G PPY(CD/K) such that CDyj = for all 



G a (K;L) :={b G G a {K) \Lb = 0}=B. 




1 i ^m. 



(ii) 
(iii) 



P\(CD/K)=K( yi ,...,y m ) 6 . 
PN{CD/Kf = K 8 . 



1 i < m. 



(ii) PPY(CD/K)=K( yi ,... 

(iii) PPN (CD /Kf =K . 



,y m )A- 
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A PPV extension for the operator (D exists and is unique up to (non-unique) differential isomorphism, thus 
justifying the notation PPY((D/K). 

Definition 1.9. The parameterized Picard-Vessiot group of D over K (or PPV group) is the set of A-field 
automorphisms of PPY ("D/K) which leave K point-wise fixed, and is denoted by Ga\^('D/K). 

The solution space of T> is the A^-vector subspace of PPY {"D/K) consisting of all y G PPV ("D/K) such that 
•Dy = 0, and we denote it by Sol(©). 

One can show that Sol(£>) is m-dimensional as a Xo-vector space, and the action of GaL^D/A') on Sol(£>) 
identifies Gal A ('D/K) with a subgroup of GL(Sol(£>)). It is shown in [3] that the image of Gal A ('D/K) under 
this monomorphism is a linear differential algebraic group. 

THEOREM 1.10 (Parameterized Galois correspondence [3, Thm. 3.5]). Let (D G K[S\. LetP be the set of A-field 
extensions LofK contained in PPY{D/K), and let S be the set of Kolchin-closed subgroups H of Ga\.&{'D/K). 

Consider the maps \|/ : 5 -> P and q> : P -> S defined by y(H) := {a G PPY(Z)/K) \ aa = a\/o G H} and 
(p(L) :=Gal A (£>/L). Then: 

(i) (p defines an inclusion-reversing bijection P f-> 5, with inverse given by \|/. 

(ii) L G P is a PPV extension of K (for some operator "D') if and only if cp(L) is a normal subgroup of 
Gal A (©/,K). In this case, o(L) = L for all o G G&\a{'D/K), and the restriction homomorphism o i-> o\l is 
a surjection Gal A (2V#) -» Gal A (©7^). 

Remark 1.11. The 8-subfield K{y u . . . ,y m ) s C PPV(<D/K) defines a PV extension of (K,5) for ?D (where 
{yi,... ,y m } is any fixed # -basis of Sol(») in PPY(Z)/K)). The 8-field PY(D/K) is stabilized by Gal A (©/^), 
and the resulting restriction homomorphism is in fact an injection 

<J ^ o|pv(©/iC) : Gal A (D/Z) Gal 8 (!D/^). 

Theorem 1.12 (Cassidy-Singer [3, Prop. 3.6 (2)]). For'D G *T[5J, Gal A (2>/^) is Zariski-dense in Galg t ( ( D/K). 

2. Algorithms for first order equations 

In this section, we give three auxiliary algorithms which will be applied in the next section to compute the PPV 
group of a second order equation. The main idea for Algorithm 1 is well-known (see [3, Example 7.1]); we 
have merely modified it in order to avoid factoring polynomials into irreducible factors. Algorithm 1 also helps 
illustrate Algorithms 2 and 3. Algorithm 2 solves an analogous problem to that solved by Algorithm 1, but in a 
slightly different context. The main idea for Algorithm 2 comes from an argument given by Singer in the course 
of the proof of [3 1 , Prop. 4. 1 (4)]. Algorithm 3 computes the PPV group corresponding to the intersection of the 
PPV fields for two first order inhomogeneous equations. 

Setting. For the rest of this paper we shall use the following notation: 

(i) Kq is a differentially closed 3 r field extension of (C(t),d t ), where d t denotes the derivation 4 on C(t). 

(ii) K := Kq{x), endowed with the structure of A := {8 V , d t } -field determined by setting 8 x x = 1, d t x = 0, and 
8 x a = for all a G Kq. 

(iii) We consider the elements of K as rational functions in x with coefficients in Kq, and we denote the degree 
of a polynomial in x and greatest common divisor of polynomials in x by deg r and gcd r , respectively. We 
adopt the convention that deg t (0) = 0. 

Consider the A-field extension F of K of given by F := K(r\), where § x Tj := pr\, d t r\ := qr\, and p,q £ C(t)(x) 
satisfy the integrability condition 5 x q = d t p. We do not exclude the possibility that r\ G K or, equivalently, 
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that F = K. We wish to compute the PPV group corresponding to the first order (inhomogeneous) differential 
equation 8 X Y = r\ over F. 

We claim this is a differential algebraic subgroup of G a (Ko). To see this, consider the associated homo- 
geneous equation [29, Ex. 1.18]: let (P := 8* — pb x . One can show that, if {1,0} is a ^o-basis of Sol(!P), 
then PPV((P/F) = F(8)a- We may assume without loss of generality that 8*8 = r|, because 8 r (8 v 6) = p8 x Q 
and 8 T 6 ^ together imply that 8*8 = ar\ for some a G K£. It follows from this that 8 v (g8 — 6) = for all 
a G GalA(^P /F). The map a i— > b a := (a6 — 6) defines an isomorphism of Gal^'P/F) onto a differential al- 
gebraic subgroup of G a (Ko). From now on, we will identify GaiA (^P/F) with its image in G a (Ko) via this 
isomorphism. 

Simply put, each algorithm will use the input data to construct a system of homogeneous linear equations, 
and then construct the output from a solution for the system which satisfies a certain minimality condition. The 
construction of these systems of equations will rely on the following Lemma. 

Lemma 2.1 (Reduction Lemma). Let p,q G K, r\ G F, and (P G K[8 X ] be defined as in the preceding discussion, 
and let ^ L G Kq [d t ] ■ The following are equivalent: 

(i) Gal A (27F)CG fl (Ab;£). 

(ii) L(r\) = 5 x f for some f G F. 

(Hi) L(r\) = (h x h + ph)r\ for some heK. 

Proof. Assume that Gal A (T/F) C G a (K ;L), or, equivalently, that L(b a ) = for all a G Gal A (^P/F). Since 

a(L(Q))=L(Q)+L(b a ) = L(Q), 

we have that L(Q) G F by the Galois correspondence. Since 8 V X(6) = L(x[), setting / := L(Q) gives the impli- 
cation (i) =>- (ii). Conversely, assume that there exists / G F such that 8 V / = L(r\). Then 8 X (L(Q) — f) = 0, so 
L(Q) = f + a for some a £ Kq. This implies that L(b a ) = for all a, and we have shown that (ii) =>• (i). 

To see that (iii) => (ii), let h G A" be as in part (iii) of the Lemma, and set / := hx\. 

Now we will show that (ii) => (iii). Since ^ £ K,we have that ^ G K for any «, so ^ G If r| G 
setting h:= J - establishes our claim. If r\ f K, we have two cases to consider: either r\ is algebraic over K or r\ 
is transcendental over K. 

If r| is algebraic over K of degree m, then {1,T|, . . . ,T| m } is a basis for F as a ^-vector space, so we may 
write / = YT=o l aiT \' w ^ tn a ' G A' in a unique way. Then 8*/ = £(8*a,- + //?a,)if and, since = y 6 A", we 
have that „ , 

— = X ( 8 -v a *' + ipai)T\'~ = 8^1 + pai , 
'M (=0 

because the set {rp 1 , 1, . . . ,r| m ~ 2 } is also a basis for F as a A"-vector space. Therefore, 8 x f = 8 x (air\), and 
setting h:=a\ establishes our claim in this case. 

Now assume that r\ is transcendental over K. Then / has a unique partial fraction decomposition: 

where a, G K and Cj ^e^ G K, because the set 

«:={T!'|^0}uU{(^|^l} 

eeK 

is a basis for £(r|) as a ^-vector space. Since the image of £f under the invertible ^-linear map given by 
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multiplication by rp 1 is still a basis, and = ^ G K, we see that 

= S^ai + pa\. 

Therefore, 8 x f = h x (a\r\), and setting h := a\ establishes the Lemma. □ 

Remark 2.2. The equivalence of (i) and (ii) in the Lemma is well-known (see [8, 31]). The fact that (iii) (ii) 
was used by Singer in [3 1]. We do not know of a proof of the implication (ii) =4> (iii) in the literature. 

The following two algorithms compute the operator L G Ko[d t ] such that Gal^('P/F) ~ G a (Ko;L). Since 
G a (Ko;L) = Sol(X), its dimension as a C-vector space is equal to ord(X), and L is the operator of smallest 
order such that the equivalent conditions of Lemma 2.1 are satisfied. Algorithm 1 applies only when r\ G K, 
while Algorithm 2 applies only when r\ ^ K. 



2.1 First case: r\€ K 

By part (ii) of Lemma 2.1, the operator L G Ko[d t ] such that GaU^/A') ~ G a (Ko;L) is the operator of smallest 
order such that Lr\ = 8 x f, for some / G K = F. We will write L and / with undetermined coefficients, and use 
this relation to obtain a system of linear equations in the coefficients of L and /. 

Algorithm 1 (Primitive of a rational). 

Input: r| G C(t)(x). 

Output: £fl := E-Lo^ G C(f)[3,] such that Ga\ A (T/K) ~ G a (K ;L^). 

Stepl: Write 



with r| f G C{t) [x] for ^ i < 2, such that 

gcd^rii,^) 



ri = no + 



Hi 
T12 ; 



1 and deg-(Tii) <deg_(Ti2) 



Step 2: Let ^ 

"~ gcd v (ri 2 ,5 r ri2)' 

and let n be the smallest integer such that d n r\ is a polynomial in x. We remark that d is the product of the 
irreducible factors of r\2, and n is the highest multiplicity of an irreducible factor of r)2- 
Step 3: Let M := deg x (d) and s := deg v (r|o). For each < N ^ M, write 



^:=£a ; 3; and /w := £ P ^ + £ L <=° 

(=0 j'=0 k=l \ 



where the a,, (3,-, and are undetermined coefficients. 
Step 4: Treating each [3 7 - and ^/ as a 8 T -constant, set 

8 x f N = L N (r\), (Hiv) 

and then multiply each side of this equation by <i n+Ar . The result will be an equality of polynomials in x, whose 
coefficients are homogeneous linear forms in the a,, py, and Equating coefficients of like-powers of x, we 
obtain a system of 
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M(n+N)+s 

homogeneous linear equations with coefficients in C(t), in the 

(N+l) + (s+2)+M(n+N-l) 

variables a,-, $j and ^ j. 

Step 5: If N = M, the system of linear equations defined by (Ha/) has a solution with not all a, being equal 
to zero. Find the smallest nonnegative integer £ ^ M such that the system of linear equations defined by (Hv), 
with N = £, has a solution with not all a, being zero. 

Step 6: Find a solution 

a,- = a t for ^ i ^ £, 

$j = bj forCKj'O + l, 

= c k i for 1 ^ k ^ n+N- 1 and < I ^ M- 1; 



for this system of homogeneous linear equations, such that each ai,bj,Ckj € C(t) and ag = 1. This condition 
determines ao,...,a,£ uniquely. 

Step 7: Set := £f =0 a;3j, and go to Output. 

Remark 2.3. Upon inspection, we see that the system of linear equations of Step 4 has coefficients in the smallest 
algebraic extension L of C(t) such that r\ € L(x). Therefore, G L[d t ]. 

Remark 2.4. In practice, it is often possible to calculate more directly, as follows: let {d\ cIm} be the set 
of poles of T|, and let n be the maximum order of r\ at any of these poles. If 

n n M 



1 = IV j +II (j[ _,u 

;=0 k=ll=l V A Ul > 

is the partial fraction decomposition of r), it is shown in [3, Example 7.1] that G a (Ko;L^) is the C-span of 
{ e u} l iL\ m ^o- Incidentally, this proves the claim made in Step 5, that the system obtained from (H,v) with 
N = M has a solution with not all a,- being 0, because 

M ^ &\m c (G a {Ko;L^)) = ord^). 
If {e~kYk=\ * s a basis of G a (Ko;L^) considered as a C-vector space, then X/ 1 is given by Y!i=Q a i^'f where 

/y(0) ... Y M\ 

e~\ d t e~\ ■ ■ ■ dfe-i 



det 



5>r (i) . 



!=0 



\ en d t e e ••• dfeej 

Algorithm 1 circumvents the need to factorize the denominator of Tj, because performing such a factorization 
could be computationally infeasible in practice. However, when such a factorization is available, it should be 
more efficient to carry out this simpler algorithm, instead of Algorithm 1 . 

Example 1. We will now apply Algorithm 1 with r\ := ^M^, to compute the PPV group corresponding to 

dY/dx = 2t/(x 2 + t). (2.1) 

In this case, d = x 2 +t, n = 1, M = 2, and s = 0. We write the operator and the rational function fy with 
undetermined coefficients, as in Step 3, with N = I: 

A:=cci| + oc and f { :=p + pix+ '-= . 

x z + 1 

Step 4 then requires us to substitute these expressions in (Ha/), with N = 1, to obtain 
3/i o , -2^ Q x-2^ A x 2 2m +2oc 1 -2mi 

-=r~ = Pi + -s- — 7-5 — — = o V-r^> — v) = £iv\)- ( 2 - 2 ) 

dx r x 2 + t (x 2 +t) 2 x 2 +t (x 2 + t) 2 y y 
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After multiplying by d 2 = [x 2 + t) 2 on both sides, we have that (2.2) holds if and only if the following polynomial 
in x, is zero 

(-Pi)* 4 + (2m + 2<xi - 2f pi +^\,\)x 2 + {2^i fl )x + (2f 2 oc - f 2 pi - f£i,i) . 
Thus we obtain a system of linear equations by setting each coefficient equal to zero. We then find the solution 
OCo = —51 0Ci = 1, Pi = Po = ^1,0 = 0, £1,1 = — 1> and we thus obtain an operator L £ C(?)[4] and a rational 
function / £ C(x,t), defined as L := ^ — ^ and / := x j^, which satisfy condition (ii) of Lemma 2,1. One can 
check that there is no solution with (Xi = 0, and therefore the PPV group of (2.1) is G a (Ko; 4 — 4). 

Proposition 2.5 (Algorithm 1 is correct). Suppose thatr\ £ K. Then, the Output of Algorithm 1 is the operator 
which defines Ga\ A (2/K) as a subgroup ofG a {K ). In other words, Gal A (P/K) ~ G a {K ;L^). 

Proof. The solution found in Step 6, for the system of linear equations of Step 4, also defines a rational function 



s+\ (n+e-i) 
j=0 k=\ 



d k 



By construction, £ Ko[d t ] and f £ K satisfy part (ii) of Lemma 2.1. Therefore, Ga\ A (!P/K) QG a (Ko;L^). 

We need to rule out the possibility that there exist L £ Ko[d t ] and / £ K satisfying part (ii) of Lemma 2.1 
with £ := ord(Z) < ord(Z T| ) = I. We proceed by contradiction: assume that such ~L and / do exist. Write 

I ,v+l T fyM-i~ i\ 

2:=£fi#, and /:=£V + £ . 

(=0 7=0 k=i \ " / 

for some T £ N, with di,bj,dk : i £ Kq. A calculation shows that d n+e L(r\) £ Kq[x\. Therefore, d n+i 8 x f £ Ko[x], 
and we may take T = n + i— 1. But this means that a,- = a,, P 7 = ^ = c^.i is a solution for the system of 
equations defined by (Ha/) with N = I, which contradicts the choice of I. □ 

2.2 Second case: r\ £ K 

By part (iii) of Lemma 2.1, the desired L £ Ko[d t ] is the operator of smallest order such that L(t[) = (b x h + ph)r\ 
for some h £ K. We will write L and h with undetermined coefficients, and use this relation to obtain a system 
of linear equations in the unknown coefficients of L and h. 

Algorithm 2 (Primitive of an exponential). 



Input: p,q £ C(t)(x) such that d,p = 5 x q. 

Output: £J>>« := ZLo^i £ Wjfit] such that Gal A (27F) ~ G a (K ;£P>' 



Step 1: Write pi ^ 

P = Pq-\ and q = q -\ , 

P2 qi 

with pk,qk G C(f)[x] for ^ k ^ 2, such that 

gcd x (/?i,p 2 ) = 1 =gcd v (<7i,<7 2 ), 
de g* (Pi)< de §.v ipt) and de g.v (? 1 ) < de &v (12 ) • 

Step 2: Let ;? 2 g 2 

gcd x (p 2 ^ 2 ,8 x (p 2 ^ 2 )) ' 
let V be the smallest integer such that cf v /? and d v q are both polynomials in x, and define 

m:= deg x (d), n:=max{deg x (po),deg x (qo),\}, and A/ := mn + n+ 1 . 

We remark that d is the product of the irreducible factors of p 2 g 2 , and V is the highest multiplicity of an irre- 
ducible factor of P2 or g 2 . 
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Step 3: Define a sequence {/?,} C K recursively: Ro := 1 and Rj := d t Ri-i +qRi-\ for i ^ 1. 
Step 4: For each nonnegative integer N ^ M, consider the inhomogeneous first order differential equation 
(with undetermined coefficients (Xo, . . . , OCjv): 

5 x Y+pY = j^a i R i . (2.3) 

Find bounds 3 Sn, G N, depending only on p, q, and N, such that the following condition holds: if ao, . . . G 
#o is any {N + l)-tuple of elements, and G K satisfies the equation obtained from (2.3) after replacing a, by 
a, for each ^ i ^ N, then must be of the form 

S N T n ( ym-l e I \ 

7=0 *=1 V / 

for some (3o, . . . , $s N Xi,o, ■ ■ -&T N ,m-\ S Ko. That is, the degree of the polynomial part of is bounded by Sn, 
and the order of any pole of % is bounded by T^, regardless of the values of ao, . . . , a/v in Kq. 
Step 5: Replace Y with this expression for in (2.3), to obtain 

N 

5 x h N + ph N = a i R i > ( I A f ) 

1=0 

where the Py and ^ / are to be treated as undetermined Sjc-constants. Multiplying by d TN+n on each side, we 
obtain an equality of polynomials in x, whose coefficients are homogeneous linear forms in the a,-, Py, and ^ /. 
Equating like-powers of x on both sides, we obtain a system of 

m(T N + n) + S N + n 



homogeneous linear equations with coefficients in C(t), in the 

(N+\) + (S N +l) + mT N 

variables a,-, Py, and ^,1- 

Step 6: If N = M, the system of linear equations defined by (1^) has a solution. Find the smallest nonnega- 
tive integer £ ^ M such that the system of linear equations defined by (I,y), with N = £, has a solution. 
Step 7: Find a solution 

a, = a t for ^ i < £, 

£>k,l = Ck,i for 1 ^ k ^ 7/v and ^ / ^ m — 1 ; 



for this system of homogeneous linear equations, such that each a^bj^c^i G C{t) and ae = 1. This condition 
determines ao,...,ag uniquely. 

Step 8: Set LP A := £f =0 a,3j, and go to Output. 

Remark 2.6. Upon inspection, we see that the system of linear equations of Step 5 has coefficients in the smallest 
algebraic extension L of C(t) such that p,q G L(x). Therefore, L p,q G L[d t ]. 

Remark 2.7. As in [31, Prop. 4.1(4)], the operator L p ^ produced by Algorithm 2 must be nonzero. Otherwise, 
the element h G C(t)(x) obtained from (2.4) by substituting [3 7 = bj and ^ = cu,i would have to satisfy the (In) 
with N = and ao = 0, which is 8 x h = ph. This implies that h = ar\ for some a G Ko, which would contradict 
the assumption that r\ ^ K. 

5 The computation of these bounds is elementary, but somewhat lengthy. We will indicate how to compute these bounds in Remark 2.9. 
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Example 2 (The Picard-Fuchs equation for the Legendre family of elliptic curves). Consider the first order 

inhomogeneous equation _ , . , . 

6 4 dY/dx = (x(x- l)(x-t)) 2 . (2.5) 

We will take r\ := (x{x — l)(x— t)) z , and apply Algorithm 2 with inputs 

d-n/dx 1/1 1 1 \ dn/dt 1 / 1 

' + r + and q: 1 



r\ 2 \x x — l x — t J " r\ 2\x — t i 

to derive the Picard-Fuchs equation for the Legendre family of elliptic curves E t :y 2 = x(x — 1) (jc— i) (cf. [24, pp. 
77-78] and [12, Example 6.9]). 

In Step 2, we compute d = x(x — 1)(jc — ?), V = 1, m = 3, n = 1, and M = 5. Now let 

be the sequence defined in Step 3 (we will have no need to compute /?; for i > 2), and consider the inhomogeneous 
equation with undetermined coefficients of Step 4, with N = 2 



37 1/1 1 1 \ „ ai / 1 \ 3a 2 / 1 



dx 2\x x—\ x — t J " ' 2 \x — t J 4 V — 2 , 
A computation shows that we may take T 2 = 1 and S2 = 1 (see Remark 2.9, or the proof of Proposition 3.5 
below). Now we write the rational function h% = Po + Pi* + 737 with undetermined coefficients po, Pi , ^0.1 as in 
Step 4, to obtain (I#) with Af = 2 as in Step 5: 

I (Po + Pi* + |f) - \ {I + ^ + (Po + M + B) = «o + f (£) + ^ ( ^) . (2.6) 

After multiplying by x(x— 1) (x — t) 2 to clear denominators, 6 and then subtracting the left-hand side from the 
right-hand side, we obtain that (2.6) holds if and only if the following polynomial in x is zero: 

(ceo + ^pi)x 4 + ({-2t - l)oo - 5CX1 + f Po - ^Pi)* 3 

+ ((t 2 + 2t)a + {\t + \)ai - la 2 - (f t + l)p - ^Pi + f^Lo)* 2 

+ (-f 2 a - \ta { + \a 2 + it 1 + §f )Po + ^ 2 Pi - (f + 2)5i, )* + (-^ 2 po + 5^1,0 )*°- 

Setting each coefficient equal to zero yields the system of linear equations defined in Step 5. One can check that 
OCo = — \, «i = 2t — 1, 0C2 = t(t — 1), Po = \{t — 1), Pi = 5, and c^o = \t(t — 1) is a solution, 7 and that indeed 

t{t - l) W + {2t - l ^-^ = Tx{{2 {t - l) + 2 X+ ^P)- 

After checking that the system of linear equations does not have a solution with a 2 = 0, we conclude that the 
PPV group corresponding to (2.5) is G a (K ; t(t- 1)^ + (2t- 3). 

Proposition 2.8 (Algorithm 2 is correct). Suppose thatr\ £ K. Then, the Output of Algorithm 2 is the operator 
which defines GaX^/F) as a subgroup ofG a (K ). tn other words, G&\ A (¥/F) ~ G a (K ;L p q ). 

Proof. Consider the sequence {/?,} C K defined in Step 3. Since 3,r| = qr\, we have that d' t r\ = /?,T| for each 
i > 0, and therefore / n \ 

i=0 \i=0 / 

6 Note that here the algorithm calls for multiplying by d 2 = (x(x— 1) (jc — f)) 2 , but in this concrete example we see directly that this is 
unnecessary. One can use square-free factorizations to do this in general (see Section 4). 

7 We have taken 0C2 = t(t — 1) instead of Ch = 1> because this is how the Picard-Fuchs equation usually appears in the literature. 



12 



Computing Galois groups of parameterized second order equations 



for each ^ N ^ M. Note also that the solution found in Step 7, for the system of polynomial equations defined 
in Step 5, also defines a rational function h G C(f)[je], obtained from the rational function (with undetermined 
coefficients) hi of (2.4) by setting = bj and ^,1 = Ck,i for each j,k,l: 



By construction, L p - q G Ko[d t ] and hi G K satisfy part (iii) of Lemma 2.1. Therefore, L p - q {b c ) = for every a G 
GaU^/i 7 ), where b a is defined as in the discussion preceding Lemma 2.1. We need to rule out the possibility 
that there exist L G Ko[d t ] and h G K satisfying part (iii) of Lemma 2. 1 with 

£:=ord(2) <ovd(L M ) = l. 

We proceed by contradiction: assume such L and h do exist. Let L := LLo^!' with d{ G Ko. By Lemma 2.1, 

h satisfies the inhomogeneous equation [ 

8 x h + ph = J £a i R i . (2.7) 

!=0 

If the bounds Si and Ti of Step 4 exist, then h must have the following form 

+ (2-8) 

j=0 k=l \ " / 

This implies that setting 5,- = a,, £y = (3/, and c^; = / for each i,j,k,l is a solution for the system of homoge- 
neous linear equations defined by (I/v) with N = £. This contradicts the choice of £ 

It remains to be shown that the bounds Si and T£ of Step 4 exist. We will compute these bounds directly, 
using a modification of the argument given in [30, Lemma 3.1], and we will see that they depend only on p, q, 
and £. Let e G be a pole ofh£ Ko(x). We claim that e must be a pole of either p or q. To see this, suppose 
on the contrary that e is a pole of h, but not of p or q. If — T is the order of h at e, with T ^ 1, then the order of 
5 x h + ph is — T — 1 (because e is not a pole of /?), while the order of Yfi=o^>^i i s nonnegative (because e is not a 
pole of <7); this gives a contradiction. Now let e be a pole of p or </, and let 

I p = — + (higher order terms) 

S ' R ' = i + ( hl g her order terms) _ * Wj 
i=0 \ x e ) h = -7f + (higher order terms). 

(x-e) T 



be the (x — <?)-adic expansions of these elements, with W2 G C{t) , and Wi,W3 6^ X . If we substitute these 
expressions in (2.7), we obtain 

(— T\\>3 \ / VV2W3 \ Wl 
7 + (higher order terms) + r^- — h (higher order terms) = r (- (higher order terms), 

(x — " / \(x — e)-'+ T2 " / (x — e) Xl 

(2.9) 

We see that either max{T + 1, T + X2} = Xi, or else the lowest order terms in the left4iand side of (2.9) must 
cancel. In order for this cancellation to take place, it is necessary that 7 + 1 = T2 + 1 and TW3 = W2W3, or, 
equivalently, X2 = 1 and T = W2- In any case, we have that T ^ max{w2,Xi — 1} (if m>2 is not an integer, we 
deduce that T ^ Xi — 1). Evidently, W2 depends only on p. Although Xi depends on the du as well as on q and I, 
one can easily find a (generic) upper bound for Xi which depends only on q and I, and not on the specific values 
of the d\\ namely, take Xi := min,{ord e (/?,■) | ^ i ^ £}. Now we may take Ti to be greater than any of the finitely 
many bounds so obtained for each pole of p and q. 

The bound Si is computed as follows: let S be an upper bound for the degree of the polynomial part of each 
Ri, for ^ i ^ I. It follows from (2.7) that maxjS'i— l,deg r (po) + Se} ^ ^> where po is defined as in Step 1. 
This concludes the proof that the bounds 7| and 5^ of Step 4 exist. □ 
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Remark 2.9. In the preceding proof, we showed how to compute the bounds Sn and Tn of Step 4 by working 
with the (jc — e)-adic expansions of p and q, for each linear factor (x — e) of d over C(t) [x]. In practice, it may not 
be feasible to compute a factorization of d into irreducible factors over C(t), let alone linear factors over C(t). 
We will now indicate how to obtain the bounds Sn and Tn without factorizing d. 

Let (ooj • • • ,fl/v) be any (N + l)-tuple of elements of Kq, and let h £ Kbea solution for the equation obtained 
from (2.3) by replacing a, with a,- for each i: 

b x Y+pY = '£a i R i . (2.10) 

(=0 



If (jc — e) is a factor of d in C(t) [x], the argument given in the proof of Proposition 2.8 shows that h can have a 
pole at e of order at most Xj — 1, where Xi is defined as in the proof of Proposition 2.8, unless p has a pole at e 
of order exactly 1. Given this a priori bound, it suffices to bound the order of a pole of h at the simple poles of p. 
We now compute the divisor of simple poles of p: 

_ = gcd x (d,p 2 ) 
gcd x (d,p 2 ,8 x p 2 )' 

This is the product of all linear factors of pi which have multiplicity one in p 2 . We will mimic the argument 
given at the end of Proposition 2.8, using powers of d\ as denominators, instead of powers of (x — e). Write 

N 

V-i Wl W2 W3 

) ajRj = — f + (irrelevant terms), p = h (irrelevant terms), and h = —~ + (irrelevant terms), 

£o d \ rf l rf i 



where W2 € C(f)[jc], wi,W3 6 Zb[je], an d deg A .(w J ) < deg t (rfi) for each 1 ^ j ^ 3. A term in £,£?i/?;, p, or /z is 
deemed "irrelevant" if its order at every linear factor of d\ is higher than — x, —1, or —T, respectively. In order 
to bound T, we proceed as before: substituting these expressions in (2.10), we obtain 



-TwT,5 x di \ I W2W3 \ wi 

r+1 h (irrelevant terms) + r+1 + (irrelevant terms) = —+ (irrelevant terms). (2.11) 



It follows that T + 1 ^ x. If this inequality is strict, then the lowest order terms in the left4iand side of (2.11) 
must cancel. For such cancellation to take place, it is necessary for (W2W3 — ^38^1) to be divisible by d\ in 
Ko[x]. This is almost what we want, except that this condition on T depends on W3, which in turn depends on the 
cij. However, since deg v (w3) < deg x (d\), we have the implication 

di\(w 3 (w2-T8 x di)) ^ gcd x (d u (w 2 -T5 x di)) / 1. 

If we now perform the Euclidean algorithm to compute this gcd,., we will obtain a series of residues. These 
residues will be polynomials in x, whose coefficients are polynomials in C(?)[r]. At least one of these residues 
must be zero in order for the gcd v to be different from 1. Therefore, T must satisfy at least one of a finite 
number of explicitly constructed polynomial equations with coefficients in C(t), and we may obtain a bound 
on the size of any integer solution for a given polynomial equation, depending only on the coefficients of the 
given polynomial. Having thus obtained the bound Tn, we obtain the bound Sn by comparing the degrees of the 
polynomial parts of the right- and left4iand sides of (2.10), after replacing Y with h and clearing denominators. 

2.3 The PPV group of an intersection of PPV fields 

In this section we give an algorithm to compute the PPV group corresponding to the intersection of two PPV 
fields of the form considered at the beginning of this section. We will apply this algorithm in §3.4 to recover the 
PPV group corresponding to a second order equation from a pair of auxiliary differential algebraic groups. 

Let r|i,r|2 £ C(t)(x), and consider the first order inhomogeneous equation 8 X Y = r\ r , for r = 1,2. If we let 
Pr '■= the discussion at the beginning of this section shows that the corresponding homogeneous equations 
are given by the operators 
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¥i:=& x -pi& x and Qz:=%-p2$x, 
and that Gal A (¥ r /K) is a differential algebraic subgroup of G a (Ko). LetE denote the PPV extension of PPV(fPi /K) 
corresponding to the operator !?2, and let 61 , 62 G E such that 8 x 6 r = T| r . Then, we have that 

K(Q r ) A = PPV(P r /K) C £ = £(6i,e 2 )A, and we may define L := PPV^i/JsT) nPPV(3 , 2 /^) C E. 

The group GaiA^r/L) is normal in G&\ A { r P r /K), since the latter is abelian. It follows from the Galois correspon- 
dence (Theorem 1.10) that L is a PPV extension of K (for some operator). If we denote the corresponding PPV 
group by A, the Galois correspondence implies that we have surjections Ji r : G&l^Tr/K) -» A, and that 

Gs1a(E/K) ~Gal A (fPi/r) XAGalA^/tf), 
where the fibered product 8 is taken with respect to the differential-algebraic homomorphisms n r . 
The following algorithm computes A, as well as the maps %\ and 712. 

Algorithm 3 (Intersection). 



Input: r|i,r|2 G C(t){x). 

Output: A differential algebraic group A 111 ' 712 , and differential algebraic homomorphisms 
7l? | Tl2 : GalA^) -» A 111 ' 112 , for r = 1,2, defined by 

A 111 ' 112 := G a (K ;L") and ^ : o ^ L' r (o) , 

such that the PPV group of PFV(&i/K) nPFV(fP 2 /K) is isomorphic to A ni,T|2 and 7I? 1 '" 12 
is the corresponding surjection of PPV groups (cf. Theorem 1.10). 



Step 1: Apply Algorithm 1 with inputs r|i and r|2, and let L\ := L"^ 1 and Lq_ := L^ 1 be the Outputs. If 
ord(Xi) ^ ord(£2), proceed to Step 2. 

Otherwise, reverse the roles of rj 1 and r(2, so that ord(iLi) ^ ord(X2), and then proceed to Step 2. 
Step 2: Write (0) (0) 

with 9 Tlf^Tlf G C(0M for ^ i 2, such that 

gcd^TlS^Tlf^^gcd^Uf), 

deg^riS^) <deg 1 .(ri ( 1 2) ) and deg^ri^) < deg Y (rif } ). 

Step 3: Let 

a := — 



gcd r (ri ( 1 2) ri^ 2) , 8 Y (ri^ri^ 

an let n G N be the smallest nonnegative integer such that d"r\ r is a polynomial in x, for each r = 1,2. We remark 

(2) (2) 

that d is the product of the irreducible factors of r\ \ and r| 2 , and n is the highest multiplicity of an irreducible 
factor of rff' or T|£\ 
Step 4: Set 

4:=ord(£,), v:=4-4, * := maxjdeg^Tif^deg^f )} , and M:=deg x (d). 
For each v ^ N ^ £\, write 

JV iV-v s+1 (n+iV-1) /rM-lt m\ 

^^Lcc-aj, 4 2) :=IM/, and .fr + I 1 



i=0 y=0 4=0 /=! 



8 Recall that this is the subgroup of elements (01,02) £ Gal^'Pi/K) x GaiA^/^ ) such that 7t; (01) = 712(02) 
9 The parenthetical superscript in r\j is an index, and does not denote an / th derivative. 
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where the a,-, |3/, Yfc> and m are undetermined coefficients. 
Step 5: Treating each yt and ^/ ,„ as a 8 v -constant, set 

£$\t\i) -tg\l\ 2 )=&xfN, (J*) 

and then multiply each side of this equation by d n+N . The result will be an equality of polynomials in x, whose 
coefficients are linear forms in the a,, (3,-, yt, and i;/ ; „,. Equating coefficients of like-powers of x, we obtain a 
system of 

J M(n+N) + s 

homogeneous linear equations with coefficients in C(t), in the 

(N +\) + (N -V + \) + {s + 2) +M(n+N - I) 

variables a,, f3/, y k , ^/, m . 

Step 6: If N = £\, the system of linear equations defined by (J#) has a solution with not all a,- being zero 
and not all P ; being zero. Find the smallest nonnegative integer CO ^ i\ such that the system of linear equations 
defined by (J/v), with N = (O, has a solution with not all a, being zero. 

Step 7: If co = l\, set L[ := L\, l! 2 := Li, and L" := 1, and then go to Output. 
Otherwise, find a solution 

a, = at for < i ^ CO, 

(3j = £» 7 for < j ^ co - v, 

Yk = Ck for ^ k ^ ^ + 1 

\im = ei.m for 1 < Z ^ n+N— 1 andO < m < M — 1; 



such that each ai,bj,Ck,ei itn G C(f) and a m = 1. This condition determines • ■ -j^ra an d • • • ,b m -v uniquely. 

Step 8: Set L[ := a,3j and := L™T V ^3/, let X" e C(f)[df] be the unique monic operator such that 10 
Xi = L"L[, and go to Output. 

Example 3. We will now apply Algorithm 3 with inputs r|i := and r\2 '■= ^r t ', note that rji = \tx\i + \. 

The first step is to apply Algorithm 1 with inputs r|i and r(2, to obtain the operators L\ := — ^ and Lq_ := 
Jf ~~ Tt ( c ^- Example 1). The polynomial d of Step 3 is given by d = x 3 + tx. We then compute the integers 
n = 1, l\ = 2, £ 2 = 1> V = 1, M = 3, s = 0, and write down the pair of operators and the rational function with 
undetermined coefficients of Step 4, with N =1: 

, (1) . ,(2) n a r £l.0 + £l, \X + Z,1,2X 2 

L l>:=aid t + a , L{>:=$ , and /1 :=Yo + Yi*+ — rrz — • 

x i + tx 

Then, we substitute these expressions in equation (J iV ) as in Step 5, to obtain 

„ 9 / x 2 +t 2 x+t \ , n x 2 +t 2 x+t d It _ d 1 „. 1 ^l.oHl.i-*+^i,2-v 2 \ , 9 ]9 n 

a l37 .x-3+tt J+«0- . t3+tt -Po-^- (Yo+Yl* + J • (2-12) 

After expanding this out and doing some simplification, we obtain 



-ajt 2 x 2 1 a o x 2 +(a t 2 +2ait)x+a o t _ 2f|3o _ y , 5i,i+2£i,2* _ (^i,o+^u*+^i.2* 2 )(3* 2 +0 



After multiplying by d 2 = (x 3 +tx) 2 on both sides, and then subtracting the left-hand side from the right-hand 
side, we obtain that (2.12) holds if and only if the following polynomial in x is zero: 

Yix 6 -a x 5 + (-f 2 a -2mi — 2rp + 2?yi -^ 2 )x 4 + (-2m -2^i,i)x 3 

+ (-t 3 a -t 2 ai -2t 2 $ +t 2 y 1 -3^ L0 + t^^ 2 )x 2 -t 2 a x-t^ lfi x°. 

10 The ring C(t)[d t ] is a left and right Euclidean domain, and there are algorithms to compute this L" (see [29, §2.1]). 



16 



Computing Galois groups of parameterized second order equations 



Setting each coefficient equal to zero yields the system of homogeneous linear equations defined in Step 5. 

|, Yi = 0, Si,o = 0, $i,i=0, $1,2 = -^ 



We verify that oco = 0, (Xi = 1, Po = \, Yl = 0, $io = 0, $u = 0, $12 = —\t is a solution, and that indeed 



j^i ~~ | r l2 = ( 2(x*+tx) ) • ^ ne sma llest value for N allowed by Step 4 is N = 1, since we assume /V ^ V = 1. 
Therefore, CO = 1, and we have L[ = 4, = I, and L" = 4 — 4. 

Our proof that Algorithm 3 gives the right answer will rely on the following well-known result. 

Theorem 2.10 (Kolchin-Ostrowski [19]). Let E be a 8 v -heid extension of K such that E 8 * = K , and let 
ei , . . . , c m , f 1 , . . . , f n G E such that ^ 6 K for each 1 ^ i ^ m and 8^f j G for each 1 ^ 7 ^ n. 

Then, there exists a nonzero polynomial G G . . . ,X m ,Yi, . . . ,Y n ] such that G(ei, . . . , e m ,fi, . . . ,f„) = if 
and only if at least one of the following holds: 

(i) There exist integers m,- G Z, not a]] zero, such that fl/li e i"'' ^ 

(ii) There exist elements Cj G K 5 *, not all zero, such that L'j=i c jfj G K. 



Proposition 2.11 (Algorithm 3 is correct). Letr\i,r\ 2 G C(t)(x), and letA^ 1 ^ 2 and 7i? 1,Tl2 be the Outputs of 
Algorithm 3. Define the operators G K[8 X ] for r = 1,2, as well as the fields E and L as above, that is 

E :=PPV(fP 2 /PPV(!P 1 /^)), and L := PFVfa/K) OPPY^/K) C E. 

Then, 71? 1 ' 112 : Ga\A('P r /K) -» A* 11 ' 112 is a surjection, and L is the PPV extension corresponding to ker(7t? 1,Tl2 ) by 
Theorem 1 . 10 applied to PPV \(B r /K), for each r = 1,2. 

Proof. Let 8 r G £ such that 8*8,- = ri r . Then, PPV{<P r /K) = K(Q r ) A , and since 5 x d?Q r = 3"rj, G if, it is enough 
to take derivatives with respect to d t only; that is, PPV ((P r /K) = K(Q r )^ t . If we define l r as in Step 4, Lemma 2.1 

impUeS ^ VPV(!P r /K) = K(Q r ) dt = K(Q n d t Q r , . . .,$-%), 

and the transcendence degree of PPV{¥ r /K) over K is £ r . 

The mapoH- (oQ r — Q r ) identifies GsA^Tr/K) with G a (Ko;Lr), for each r = 1,2. By Theorem 1.10, there 
is a unique 3 r algebraic subgroup H r ^ Gdl A { r P r /K) such that L is the fixed field of H r ; that is, L is precisely the 
set of elements a G PPV such that at? = a for every a G //,-. On the other hand, if L, L' G Ko[d t ] are such 
that G a (£o;L') Q G Q (^ ;^), then there exists X" G Ko[3 f ] such that L = £,"£,', and G a (£b;£)/G a (Kb;£') ~ 
G a (K ;L") (see [29, p.45]). If L' r G XbPt] is a right-hand factor of L r , the fixed field of G a (K ;L' r ) is precisely 
K(L' r (Q r ))A- Therefore, there is a right-hand factor L' r of L r such that L = K(L' r {Q r )) A . 

We remark that every element a G L\/T defines an algebraic relation, as follows: since a G PPV{¥\/K), there 
is a rational function <I> G K(Y\, . . . ,Y^) such that a = <&(Qi, . . . ,d c t l ~ i 9 t ). Similarly, there is a rational function 
*P G K(Y\ ,Yg 2 ) such that a = *P(02, • • • ,3 ( 2_ 8 2 ). Setting these two expressions for a equal to each other and 
clearing denominators, we obtain that L 7^ K if and only if there is a nontrivial ^-algebraic relation amongst the 
elements 81, 8 2 , 3 f 8i, d t Q 2 , 3f 1-1 8i, 3f 2_1 8 2 €E. 

By Definition 1 .8, we have that E s * = Kq = K^ x , and we may apply Theorem 2.10 with 

fa/8i forO^ j'^i-1, and 
lar 4 8 2 for^ 



h 



Therefore, if L 7^ K, we have that there exist elements oq, . . .a^-ijbo, ■ ■ ■ ,^4-1 G Kq, not all zero, such that 

£' 2 -l 

£ a$9i - £ V/82 G if. (2.13) 

i=0 ;=0 

There is a nonnegative integer 0) G N which is minimal with respect to the property that (2.13) holds for some 
dj,bj G Kq, with a,- = for all i > cb, and aa / 0, and if we impose the further condition that a& = 1, this deter- 
mines oq, .. . ,a&,bo, •• • ,bg 2 -\ uniquely. The minimality of ct> implies that L is A-generated over K by £, 0^61 
(and therefore also by X^5/3/0 2 ). 
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We claim that go = GO, where go is defined as in Step 6. Note that this holds if and only if a; = a,, and bj = bj, 
where a^bj G C(t) are the elements found in Step 7. For definiteness, let us set 

£ai3|9i - £ bjdiQ 2 =feK, (2.14) 

r=0 7=0 

where a® = 1. If we define Zj := L£Loa;3}, := £;=0 fy^f > an ^ 

Jt r :Gal A (iP,/^)— ^G 8 (Xb) 

a i ^z;(ae,-e r ) 

for r = 1,2, we have that 7li (GalA(fPi/^T)) = %2(Gdl&(!p2,/K)'), as differential algebraic subgroups of G a {Ko). 
If we let L" G Ko[d t ] be the unique monic operator such that 71 1 (GalA^i/^T)) = G a (^o;Z"), it follows that 
L"L[ = L\, and therefore t\ — GO = ord(Z") = ^2 — ord(Z^), or in other words 

orA{U 2 ) = &-{£ l -l 2 ). 

Therefore, bj = for all j > GO — V, and Sa-v 7^ 0, where v := l\ — £2, as defined in Step 4. This implies that, 
if we apply S v on both sides of the equality in (2.14), we obtain that do, . . . ,aa,,bo . . . ,ba- v , together with the 
coefficients of a base-<i expansion of /, define a solution for the system of equations defined by (J a-), with tb = N. 
Therefore, go = GO and we obtain our result. □ 

3. Algorithms for second order equations 

Let K = K{)(x) be the A-field defined at the beginning of §2. In this section, we will describe an algorithm to 
compute the PPV group corresponding to a parameterized differential equation of the form 

5 2 x Y + n 8 x Y + r 2 Y = 0, (3.1) 

with r,- G C(x,t) for ^ i ^ 2. After performing a change of variables (see [29, Ex. 1.35 (5)] or [22, p. 5]), we 
may put this equation in the following form 

5 2 x Y-rY = 0, (3.2) 

where r := \r\ + \h x r\ - r 2 . 

The advantage of performing this change of variables is that, if 'D := 8^ — r, then Ge&$ x ('D/K) is an algebraic 
subgroup of SL2(^o) (and therefore Gs&a(2)/K) is a differential algebraic subgroup of SL 2 (Kq)), so that there 
are less possible candidates for (parameterized) Picard-Vessiot groups to consider. For the time being, we shall 
limit our attention to equations of the form (3.2). In §3.4, we will show how to express the PPV group of the 
original equation in terms of the PPV group of (3.2). 

3.1 Kovacic's algorithm 

In [22], Kovacic describes an algorithm which: (i) finds all Liouvillian solutions of (3.2), whenever such solu- 
tions exist; and (ii) computes Galg (2)/K) from the data of these solutions (or their nonexistence). Kovacic's 
algorithm proceeds by first determining whether the Riccati equation [29, §4.1] for the operator "D, which is 

giVenby R„(U):=8 x U + U*-r = 0, 

has a solution u G K. If such a u exists, a computation shows that D factors over K{u) [8 V ] into first order factors, 
as <D = (8 V + u) o (S x - u) [29, Ex. 1.36 (lb)]. Therefore, Ga\{p/K{u)) is reducible [7, Prop. 2.1]. 

Definition 3.1. We denote the upper triangular and infinite dihedral subgroups of SL 2 (Kq) by UT(^o) and 
Doo(A'o), respectively. They are defined as 



UT(£n 



aeK*,beK } and D_(*b) := { ( Q fl -iJ>l_ a -l 



a \ / a 



a G K£ 
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For a proof of a significantly sharper version of the following result, see [29, §4.3.4] or Kovacic's original 
paper [22, §1.2]. We have only included those parts which are relevant for our purposes. 

Theorem 3.2 (Cases of Kovacic's algorithm). Kovacic's algorithm falls into four main cases: 

(i) IfRqj{u) = for some u G K, then Gal§ r (£>/K) is conjugate to a subgroup of UT(Ko). 

(ii) IfR^(u) = for some u G F\K, with F a quadratic extension of K, then Gal^('D/K) is conjugate to a 
subgroup ofDoo(Ko). 

(Hi) If Cases (1 ) and (2) do not hold, and Rd(u) = for some u G K, then Gsl\^ x (D/K) is Unite, 
(iv) If Cases (1), (2) and (3) do not hold, then Gal &x (0/K) = SL 2 (K ). 

Theorem 3.3 (Dreyfus [8, Prop. 9]). Suppose that Gal^T) / K) = SL 2 (K ). Then, Gal A (£>/^) is conjugate to 
SL,2(C) if and only if the following inhomogeneous equation has a solution in K: 



The following result is a corollary of Theorem 1.12. 
Theorem 3.4. Gal 8v (£>/K) is finite if and only if Gal A (CD/K) is finite. In this case, Gal A (£>/K) = Gal Sv (£>/K). 

One can show that (3.3) has a solution in K if and only if it has a solution in C(t)(x), and there is an algorithm 
to decide whether or not this holds [29, Prop. 4.1 and Rmk. 4.5 (2)]. Therefore, Theorems 1.5, 3.3 and 3.4 
completely solve the problem of determining Gal A (£>/K) if Gal 8j (©/K) is finite or if Gal 8t (?D/K) = SL 2 (K ). 

We shall now restrict our attention to cases (1) and (2) of Kovacic's algorithm. From now on, we will 
assume that Ga\ A ('D/K) is infinite, unless explicitly stated otherwise. Our strategy is similar to the one followed 
in [8], and the arguments given there were of great help to us in developing ours. Since we do not have an 
interpretation for the solutions for our differential equations in terms of path integrals, we cannot directly apply 
these results in our context. We will give differential-algebraic versions of the arguments presented in [8, §4.1, 
§4.2], in order to prove that our Algorithms 4 and 5 are correct. In cases (1) and (2) of Kovacic's algorithm, we 
use the solutions for (3.2) to define semi-invariants of the linear differential algebraic group Gal A (©/^), in the 
sense of [1, II, §7]. In the non-parameterized setting, this approach has led to generalizations and simplifications 
of Kovacic's algorithm which work over more general fields; 1 1 see for example [15, 34]. See also [14] for an 
efficient algorithm to compute first order right-hand factors of f D. 

3.2 Upper triangular case 

Suppose that Rd(U) = has a solution u G K. One can find such a u G C(t)(x) explicitly (see [14,34] and [29, 
Prop. 4.9]). 

Algorithm 4 (Upper triangular). 



& X Y -4rb x Y -2(5 x r)Y = -23, r. 



(3.3) 



Input: u G C(t)(x). 

Output: A differential algebraic group \JT(Kq;A" ,B u ), defined by 




and such that Gal A ('D/K) ~ UT(^ ;A",B"). 



Step 1: Apply Kovacic's algorithm to compute Galg ^©/^). If it is found that 




n 



Kovacic's original algorithm assumes that the ground field K is C(x). 
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for some ngN, set A" := fi n (recall fi„ denotes the group of n th roots of unity in C x ) and go to Step 4. 
Otherwise, go to Step 2. 

Step 2: Apply Algorithm 1 with input r\ := d t u, let X 11 be its Output, and set 

A u :=G m (K ;L^). 

Step 3: If / 1, let 

B" :-- 



iiGahX'D/K)~G rn (Ko)\ 
G a (Kq ) otherwise , 



and go to Output. 
If £ n = 1, proceed to Step 4. 

Step 4: Find v G C(t)(x) such that 8 x v = d t u. Such a v exists whenever A" C G m (C). 

Step 5: Apply Algorithm 2 with inputs p := — 2u and g := — 2v, and let L p,q be its Output. Set 

B u :=G a (K ;L p ' q ), 

and go to Output. 



Proposition 3.5 (Algorithm 4 is correct). Letu G C{t)(x) such thatR^u) = 0. Then, 

Gal A (2V*r) ~UT(Ko;A M ,fl M ). 

Proof (cf. [8, §2.2]). There is a basis {ji,y 2 } of Sol(©) C PPY(2)/K) such that 8^1 = uyi and 8 JC (^) = 
Vj -2 [22, §1.3]. This choice of basis identifies GalA((D/K) w i tn a subgroup of UT(Ao), as follows: for a G 
Gal A (2)/^), define a G := ^1, and b c := ^ - A computation shows that a a G K£,b G G ATq, and 

a c 

is the embedding Gal&(CD/K) <^-> SL^A^o) determined by the choice of basis {^1,^2} of Sol(©). 

Moreover, the map % : Gal^CD/K) — > G m (Ko) : o \— > a c is a differential-algebraic homomorphism (so its 
kernel is Kolchin-closed). The PPV extension corresponding to ker(%) by Theorem 1.10 is F := K(yi)&, which 
is a PPV extension of K for the operator 8 T — u. The induced character % : GaL^S* — u) /K) — > G m (Ko) is now 
injective, and its image is a differential algebraic subgroup A of G m (Ko). If a G ker(%), then b a = d(^) ~~ so 
the map a 1— >• & c identifies GaiA(2V-F) with a differential algebraic subgroup B of G u (A"o), an d we nave that 



a 
a" 1 



a eAi G £ 



Gal A (2V*) si \JT{K ;A,B) = 

We will show that A=A U and B = B". 

First we calculate A. If Gal^CD/K) is infinite, Theorem 1.12 implies that A = fi„ if and only if 

a 



GzUV/K) ~ 



a 



-1 



a *E p. n , b £ Kq 



If A is infinite, a computation shows that <j(^ i ) = + and 8 J(r (^ i -) = Applying Proposition 2.5 with 

ri = d,u and 8 = we deduce that A = A". 
1 yi 

In order to calculate the operator defining B, we have to consider the possibilities A C G m (C) and A ^ G m (C) 
separately. The following result is proved in [32] under different hypotheses, but the proof works in our setting 
as well. 

Lemma 3.6. LetB be a nontrivial proper differential algebraic subgroup ofG a (Ko), and let a G G,„(Kq). Then, 
B = aB := {ab \ b G B} if and only if d,a = 0. 



20 



Computing Galois groups of parameterized second order equations 



Proof of Lemma. Let L G Ko[d t ] be the unique monic operator defining B. We have that I/O and ord(X) ^ 1. 
If a G G m {Ko), the monic operator defining aB is L' := cT 1 o Lo a. If d t a / 0, then L ^ L', so B ^ aB. □ 



Case 1: A(f- G m (C). The group A acts on B via 



c? 

a' 1 



1 ft 
1 



1 



1 a 1 



a) V0 1 

By Lemma 3.6, in this case B is either or & a (Ko), (because a G C 44> a 2 G C). By Theorem 1.12, B = if and 
only if Gals^/tf) ~ G m (K ), and we have that B = B u in this case. 

Case 2: A C G m (C). In this case, it is shown in [3] that there exists v £ K such that d t y\ = vy\ and 8 x v = d t u. 
Note that this implies that 8 x (yJ" ) = — 2w;y^ 2 and 3f (y] - ) = —2vy^ 2 . Since 

PPV(©/F)=F(y 2 ) A = F<|> A 
and 8 X (^) = y^ 2 , we apply Proposition 2.8 with r\ = y\ 2 and 6 = ^ to deduce that B = B" in this case also. □ 
3.3 Infinite dihedral case 

Suppose that R^U) = does not have a solution in K, but /?©(m) = for some u in a quadratic extension of A'. 
In this case, one can find (j) G C(f)(jc) (see [30,34], [29, Prop. 4.24]) such that the minimal polynomial of u over 
^is: 



Algorithm 5 (Infinite dihedral). 



Input: $ G C(?)(jc). 

Output: A differential algebraic group D^A^A^), defined by 



Doc(^ ;A^) := 



a 
a-' 1 ' 



a 
1 



-a 



aGA* 



and such that Gal A {p/E) ~ Doc(£ ;A*). 



Stepl: Set 



Step 2: Apply Algorithm 2 with inputs p := ^ and g := and let L p ' q be the Output. 
Step 3: Set A* := G m (^ ;£ M ), and go to Output. 



Proposition 3.7 (Algorithm 5 is correct). Letu G C(x,t) andty G C(f)(x), such thatR^u) = and the minimal 
polynomial of u over K is given by 

U 2 +$U+ Q8^ +^ 2 -r). 

Then, Gal A (2V^) ~ D„(I ;A*). 

Proof, (cf. [8, §2.3]). If v is the other root of the minimal polynomial of u, one can show that there is a basis 
{ji;^} of Sol(2?) such that h x y\ = uy\ and 8 v j2 = vj2- Since K(u) C PPV(©/A'), Theorem 1.10 gives a sur- 
jection GalA(©/^) -» Ga\{K(u) / K) . After possibly replacing j2 by aj2 for some a G ^ , one can show that 
t(yi) := ^2> x (3 ; 2) := — yi> := v defines a lift x G Gal^CD/K) of the unique non-trivial ^-automorphism of 
K{u) [22, §4.3]. For a G Gal A {CD/K) we define 71(a) := if a(u) = u, and 71(a) := 1 if a(w) = v. If we define 

a<; := 
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a computation shows that a a G Kq. One can show that the embedding of GoIa('D/K) in SL2(-fiTo)> determined 
by the choice of basis {yi,y 2 } of Sol(©), identifies GoIa(CD/K) with a subgroup of D^Kq), as follows: 



OH *' 

The group GaiA^/^) has two connected components, and the normal subgroup Gal&(CD/K(u)) is the 
connected component of the identity. One can show that o(y\y 2 ) = ±yiy2 for every a G Gsl\a('D/K). Therefore, 

PPV (£>/*(«)) =K(u)(y l ) A , 

and Gal A (f /K{u)) is the PPV group for the operator 8 A — u over For a G Gal A (<D / K(u)) , we have that 
a a = ?p-. If we let A C G m (^o) be the (isomorphic) image of GoIa('D/K(u)) under the differential-algebraic 
homomorphism o i— )• a a , we have that GaiA^/A') ~ D C x 3 (^o^) (defined as in the Output of Algorithm 5). We 
will show that A = A^. 

As in the upper triangular case, a calculation shows that a(^) = ^ i + ^f a an d8 A: (^ 1 ) = Note that we 
cannot apply Algorithm 2 to compute the PPV group corresponding to the inhomogeneous equation 5 X Y = d t u, 
because the conditions %^ EK and %^ G K are not satisfied. However, one can show that (b = Y : - Q - , for some 
<?i G C(f) and c, G (see [22, Step 3, §4.1]), and a computation shows that go := G satisfies 8*00 = 

By making the substitution g := / + ^co, we see that, for any L G ^o[3 r ], 

there exists / G K(u) such that L(d t u) = 8 x f <^=^ there exists g G K(u) such that il (d,u + i3 f <|)) = 8*g. 

Now we may apply Proposition 2.8 with r\ = d t u + 2 d t § and 6 = ^ + ^(0, because ^ and ^ are in C(f)(;e). 
This concludes the proof that A = A^. □ 

3.4 Recovering the original group 

At the beginning of this section, we performed a change of variables on (3.1) to put it in the form (3.2). We will 
now indicate how to compute the PPV group corresponding to (3.1). The corresponding problem in classical 
Picard-Vessiot theory could be regarded as an exercise; indeed, Kovacic [22, §1.1], as well as numerous other 
authors, consider the problem solved after computing the group corresponding to (3.2) only. In the parameterized 
situation there is a new subtlety, stemming from the fact that G m (Ko) has many infinite differential algebraic 
subgroups (cf. Theorem 1.4); while any proper algebraic subgroup of G m (Ko) must be finite cyclic. In general, 
the differential-algebraic relations amongst the solutions for (3.2) and the solution for the change-of- variables 
operator (i.e., the operator E denned below), form an obstruction to expressing the PPV group of (3.1) as an 
almost direct product (that is, a finite-index quotient of the direct product) of the PPV groups corresponding 
to these two operators. This is drastically different from the classical situation, where the PV-group of (3. 1) is 
always an almost-direct product of the PV-group of (3.2) and the change-of-variables group. 
To fix notation, consider the differential operators 

'D:=8 2 x -r, E:=b x + \r u and <J := l\ + n 5 x + r 2 , 

where r := \r\ + \h x r x - r 2 . Let E := PPV (£ /PPV (£>//<:)), and let w G E x such that Ew = 0. In other words, 
K(w) A = PPV(E/K) C E. Let y u y 2 G PPV (CD/K) be a ^ -basis of Sol(£>). If we let Zi := wy t G E, one can 
show that fzi = 0, and since zi and z, 2 are A^-linearly independent, they form a basis of Sol(^F) in E. Therefore, 

K{zi,z 2 )a = PPV(f/K) CE = K( yi ,y 2 ,w)A. 

Moreover, a computation shows that Z\5 x z 2 — z 2 5 x z\ = aw 2 for some a G Kq [29, Ex. 1.14 (5)], so either 
E = G&\a{!F /K), or else E is a quadratic extension of PPV(^F /K), generated by w. For later use, we define 
V := [E : PPY(J r /K)], the algebraic degree of this field extension. One can show that v = 1 if and only if 
Gal A {E/K) is finite of odd order. Now consider the A-subfield of £ defined by L := PPV ( t D/K) C\PPV{E/K). 
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Since Gal A ('E/A') is a differential algebraic subgroup of G m (Ko), it is abelian, and the Kolchin-closed subgroup 
corresponding to L, as a A-subfield of PPV(£/fir), is normal. Therefore, L is a PPV extension of K (for some 
operator). If we let A : = G&\a(L/K), we have a lattice of PPV extensions over K: 

E (3.4) 



PFV((D/K) 



Gal A (©/Jf) 




PPV (!F/K) 



Gal A (f/K) 



The Galois correspondence (Theorem 1.10) implies that 

Gal A (E/K) ~Gal A (-D/K) x A Gal A (E/K), 

where the fibered product is taken with respect to the surjections cp : Gd\&{'D/K) -» A and \|/ : Gal A ('E/A') 
induced from the lattice (3.4). It only remains to compute the middle term in the exact sequence 

1 — ►jiy — > Gal A (!D/^) x A Gal A (£/tf) — > Gal A (^/£) — ► 1. 



A, 



(3.5) 

The computation of Gal A (£/,£) is analogous to the computation of the diagonal part of G&\&('D/K) in the 
upper triangular case (cf. the proof of Proposition 3.5): we first compute the PV group Gal§ t ( "£/,£). If it is finite, 
then Gal A (£/A') = Gal§ x (E/K). If G&l^J^E/K) is infinite, one shows that a \- > ^ =: e a induces an isomorphism 
of Gal A ("E/K) onto a differential algebraic subgroup of G m (K ), and we have that c(^) = ^ + ^ and 
S T (^) = — \d t r\. An application of Algorithm 1, with input r\ := — \d t r\, will give G C(?)[3 ; ] such that 
Gal A ('E/A') ~ G m (Ko;L^). We will now compute A, as well as the maps cp and \|/. There are two cases to 
consider, depending on whether Gal A ('E/A') is finite or infinite. We will apply the following result in the case 
that Gal A (E/K) is finite. 

Lemma 3.8. The groups Galg r (!2) / K) and GdX^'D/K) have the same finite quotients. 

Proof. If G is a linear differential algebraic group, we denote the connected component of the identity of G by 
G°.If y:G-»H is a 3 ? -algebraic surjection onto a finite group H, then (p factors through G — » G / G" [1, p. 906]. 
Theorem 1. 12 implies that Ga\&('D/K) /Gal A (©/^)° ~ G&\CD/K) /Gal 5A (©/^)°. □ 

Case 1: Suppose that Galg r (£/A') = jd m for some m G N, so that A is also finite, the order of A divides m, 
and w m 6 K [20]. We shall list all the finite cyclic quotients of the algebraic subgroups of S\^2(Kq), in increasing 
order of the largest finite cyclic group which they admit as a quotient, and examine the possibilities for A in each 
case. 

The groups G a (Ko), G m (A"o), \JT(Kq), an d SL^^o) are connected, and the alternating group A5 is simple, 
so none of these groups admit a nontrivial finite cyclic quotient. Therefore, if Galg^©/^) is isomorphic to any 
of these groups, then A must be trivial, and we have that 

Gal A (jF/tf) ~ n mh ■ Gz\ A (<D/K) C GL^o), 

where fi m / v -Gal A ( ( D/K) := {£-c | C, € ^ m / v ,a G Gal A (©/^)} and n m i y C C x denotes the (m/v) th roots of unity. 

Each of the following groups admits precisely one nontrivial finite cyclic quotient, of order 2: the dihedral 
groups T>2n, the infinite dihedral group Dco(A'o), and the symmetric group 54. Therefore, both cp and \|/ are deter- 
mined by A in any of these cases. If m is odd, we conclude that A is trivial, and Gdl^J- /K) ~ ^ m ■ G&\&('D/K). 
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If m is even, one checks whether the (unique) quadratic extension of K which is contained in PV(CD/K) is 
(algebraically) isomorphic to K(w^) using finite Galois theory. This holds if and only if A = ^2- 

Similarly, the group A4 has a unique cyclic quotient, of order 3. If m is not divisible by 3, then A = {1} and 
we are done. Otherwise, one needs to decide whether the unique cubic extension of K which is contained in 
PV(CD/K) is isomorphic to K(w^); this will hold precisely when A = ^3. 

Lastly, for n£N, the groups of the form 



a b 
a' 1 



a G ju n ,b G B 



whether B = {0} or B = G a (Ko), each have precisely one quotient of order /, for each / dividing n, given by 
(0 cr l ) ^ a>lr ^ or some e ^ ^ relatively prime to To compute A, we proceed as before: one has to find the 
largest I G N such that: (i) I divides gcd(m,«); and (») the unique finite cyclic extensions of K of degree /, which 
are contained respectively in PV('E/K) and PV(CD/K), are algebraically isomorphic. In this case, we have that 
A = fjLi. If there is no such / G N, then A = {1}. 

Case 2: Now suppose that G&\&( ( E/K) is infinite, and therefore Kolchin-connected. In this case, A must 
either be trivial or infinite. The only differential algebraic subgroups of SLv(Kq) which admit an infinite abelian 
quotient are the groups \JT(Kq;A,B) (defined as in the Output of Algorithm 4). Therefore, A = {1} whenever 
Gsl\a('D/K) is not of this form. The following algorithm computes A and the maps cp and \|/. For simplicity, 
Algorithm 6 assumes that Gal^D/K) * s upper-triangular, and that Gsl^('E/K) is infinite and has already been 
computed. 

ALGORITHM 6 (Recovering the original group in the upper triangular case). 



Inputs: T) = 8 2 — r and £ = 8 r + \r\, with r, r\ G C(x,t). 

Output:A linear differential algebraic group A® ' E , and differential-algebraic homomorphisms 



(f>^ : Gal A (!Z)/^) -» A ' and v? v < L : Gal A CE/K) -» A 

such that A 2 ^ is the PPV group of PPY(D/K) n PPV (£/£), and (p®- 2 and are the 

corresponding diferential-algebraic surjections of PPV groups. 



Step 1: Compute a solution u G C(t)(x) to the Riccati equation Rd(U) = 0. Apply Algorithm 4 with input 
u, and let Gal A ('D/K) ~ \JT(K ;A U ,B U ) be the Output. 

Step 2: If A" C ^ 2 , compute a solution y 1 G C(x, t) to the differential equation b x Y = uY . Apply Algorithm 3 
with inputs r(i :=y^ 2 ,ri2 := -\<d t n, and let A 11 ' 112 , 7t 1 | ll!Tl2 , and^'' 112 be the Outputs. Set 



and then go to Output. 
If A u <£ /J2, proceed to Step 3. 

Step 3: If A" = fi n for some n > 3, set A® 2 - := {1}, let (p® 2 -, y®'^ be the trivial maps, and go to Output. 
If A" is infinite, proceed to Step 4. 

Step 4: If there exist integers 12 mi,rri2 G Z, not both zero, and / G C(t)(x), such that 

f ' 

there is a unique pair of integers m\ and m2 satisfying (3.6) with m\ positive and as small as possible. Set 



9 



and 



• 71 



.TI1.TI2 



m 2 

m\u — —r\ 
2 



(3.6) 



a ©,e . = | a m, J a G A »} = | a »2 J a G Gal A (£/£)}, 



and 



y 25 ^ :a^o m2 ; 



12 



See [7, Prop. 2.4] for an algorithm to compute these integers, if they exist. 
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and go to Output. 
Otherwise, proceed to Step 5. 

Step 5: Apply Algorithm 3 with inputs r|i := d t u, r\ 2 := -\d t r\, and let A 11 ' 712 , tt^ 1 ,Tl2 , and %\ ur[1 be the 
Outputs. Set 

A 2>, 2 :=a tu,H2 j y<D,<E :ai ^rf\ito ; and ^.^^ ^ 

and then go to Output. 

Our proof that Algorithm 6 gives the right answer will rely on the following result. 

Lemma 3.9. Let E := 8 A - g for some g €K, and let ^ p G PP\(E/K) be a solution (i.e., 5 x p = gp ). Let I 
be the order of the differential operator L G K [d t ] such that Gal A (E/K) ~ G m (K ;L). Then, PPN(E/K) has 
transcendence degree i+ \ over K, and it is algebraically generated by p and ^ , d t ( ^ ) , . . . , dj ~~ 1 ( ^ ) . 

Proof. Since 8 A p = gp, we have that in fact K(p)a = PPV(E/K) is generated over K by {8"p}, for n G Z^o- 
That is, PPV(E/K) = K(p) dr A computation shows that 

^(f)-f E K( P ,d t p,...,dr l p), 

and therefore K{p)^ = ^(p)(^)g ■ One can show that 8^9"(^) = 9" +1 g G ^ for every n ^ 0, and our result 
now follows from Theorem 2.10, since L is the operator of smallest order such that L(^£) E ^ (cf. Lemma 2.1 
and the proof of Proposition 3.5). □ 

Proposition 3.10 (Algorithm 6 is correct). Suppose that Gal A (£> / 'K) is isomorphic to UT(K ;A,B), and that 
Gsl\a(E/K) is infinite. Define the fields L and E as above, that is 

E:=PPV{<E/PPN{<D/K)), and L := PPV (CD / K) (~)PPV (E / K) C E. 

Then, the Output of Algorithm 6 is correct. That is, A ~ A 232 -, cp = cp 25 ' 2 -, and \\f = v]/ 25 ' 2 -. 



Proof. Let yi,y 2 G PPN (CD /K) be the solutions computed in Proposition 3.5, that is, y\ and y 2 satisfy the rela 



tions 8 x yi = uy\ and 8 Y ( — ) = y t 1 . We will divide the proof in two cases, depending on whether A is finite or 



infinite. 

Case 1: If A = \i n for some n G N, the Galois correspondence implies that y" G K, and we showed in the 

proof of Proposition 3.5 that PPY(CD/K) = K(y 1 )(^) A . Suppose thatB = G a (K ;Li), where ^ L\ G C(*)[3 t ], 

and let l\ := ord(£.i). If we let 6i := ^, the argument given at the beginning of the proof of Proposition 2.11 

shows that in fact , % , „ . , „ 

FPV(0/K) = K( yi , 6i , 3,6i , . . . , af 1 ~ 1 Gi ). 

On the other hand, if we let L2 G C(t)[d t ] such that Gal A (E/K) ~ G m (#o; la), and 8 2 := we may apply 
Lemma 3.9 to PPV(£/^) to conclude that 

PPV(£/*T) = K(w,Q 2 ,d t Q 2 , ■ ■ .,d £ r%), 
where £ 2 := ord(£2). Theorem 2. 10 shows that 

K( yi )nK(e s ,d t e s , . . .,$-%) =k = K(w)nK(e s ,d t e s , . . .,%-%), 

for each s = 1,2. Moreover, since y\ is algebraic and w is transcendental (over K), another application of Theo- 
rem 2. 10 gives that K(yi)DK(w) = K. Therefore, if we let p x := pi := ~^ L . and fP, := 8^ - p s S x for 
5=1,2 (we remark that these are the homogeneous operators corresponding to the first order inhomogeneous 
equations 8 V F = y\ 2 and 8 A F = — \d t r\, respectively), we have that 

PPV(CP S /K) = K(Q s ) dt C E, and L = PPN(<P l /K)r\PPN('P 2 /K) C E. 
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If n ^ 2, then y 1 G K, in which case Proposition 2.11 establishes our contention. We claim that if n ^ 3, 
then L = ^. To see this, we proceed by contradiction: suppose that L ^ K. By Theorem 2.10, the elements 
0i,02,3f0i,3;02, • • • ,3, 1_1 0i,3f 2_1 02 G ^ must satisfy arelation of the form 

e-i-i h-\ 

Y d a0 t Q l -Y d b j d{Q 2 =feK, (3.7) 

i=0 j=0 

with ao, ... ,ae ,-i,&o> • • • G not ai l zero. Applying 8^ on both sides of the equality in (3.7) gives a contra- 

diction, since £,a,3jrii = (L ; #,^i)) ; i 4- where the sequence {/?,} C 7? is defined as in Step 3 of Algorithm 2. 

Case 2: IfAis infinite, thenyi is transcendental over K. Let L3 €C(t)[d t ] such that A = G m (Ko;L}), and let 
6 3 := ^l. If we apply Lemma 3.9 to the field K{yi) A = PPV((8. V - u) /K) with p = y\, we obtain 

K( yi ) A = K( yi ,Q 3 ,d t Q3,.-.,dt 1 ^), 

where £3 := ord(Zg). We claim that L = K(y\) A D K(w) A C E (in other words, we may ignore the second 
solution y 2 G FPW ((D/K) in this case). To see this, let 0i := ^ and 02 := ^ as above. Since every ele- 
ment of (K(yi ,yz)A H ^(w) a)\^(vi ) a defines a nontrivial algebraic relation over K(y\)& amongst the elements 
01, 62,3,61,3,62, • • • ,3f 1_1 0i,3r _1 2 (cf. the proof of Proposition 2.1 1), and since 8^(3j0i),8^(3/0 2 ) 6 K( yi ) A , 
Theorem 2. 10 implies that, if K(w)a n ^(vi , V2)a ^(vi)a, then there is a nontrivial linear combination 

£ fli 3j0i - £ V 7 '02 =/ e ^(vi) A (3.8) 
i=o 7=0 

with G /To. Define the sequence {/?,} C A' by 7?o := 1, := d t Ri-\ — IQ^Ri-i for / ^ 1 as in Step 3 of 
Algorithm 2. Applying 8 Y on both sides of the equality in (3.8), we obtain 

£ br 2 + ^ I M/ n = (3.9) 

V S=0 / 2 7=0 

There is a unique partial fraction decomposition: / = Y*k c ky\ +Lm / 1 ~ y with Q G ^(03, • . . ,3f 3 ^3), and 
d/ j/n , e\ algebraic over this field. Therefore, if we let eo := 0, (3.9) implies that 

8v<imj(vi - e m ) - ld m j(uyi - 8 x e m ) 



8,/ = £(8,c, + ^)y* + £ 

k mj e mJ 

- lud m j ld m i (5 x e m — ue„ 

m,l 



1(8^+^ + 1 (yi _ gmy , 



-2 



= 8 v c + (8 x <i ,2 - 2ud , 2 )y 

This implies that $ x (Y,jbjd J t Q 2 + co) = 0, whence £;£/3;?02 = g — Co for some g G Xq! or > m other words, 
E ; fcy3/02 G A"(yi) A . We have shown that if Ly^-3/02 G A"(yi,y 2 )A, then E/fc/3/02 G ^(vi)a- This establishes 
our claim that L = K{y\) A C\K{w) A . 



-lfl\_r_w...\ n wa 3 a r^-li 



Since G and 8^(3j0 2 ),8 x (3/0 3 ) G Theorem 2.10 implies that 



K(y l )nK(Q s ,d t Q s ,...,d t r i Q s )=K = K(w)nK(Q s ,d t Q s ,...,d l r Q s ), 
for each s = 2,3. Therefore, if L / K, precisely one of the following statements is true: 

(i) L = K( yi )nK(w). 

(ii) L = K(Q 2 ,d t Q 2> . ■ ■ ,3? 2_1 e 2 ) nyr(0 3 ,a,03, . . .,#-%). 

By Theorem 2. 10, (i) holds if and only if there exist integers m\ and m 2 , not both zero, such that y™ l w m G K, 
and this holds if and only if the condition of Step 4 holds. In case (ii), Proposition 2.11 establishes our result. □ 
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4. Concluding remarks 

4.1 Consequences 

We mention two consequences of the algorithms. 

Theorem 4. 1 . Every unimodular PPV group of a second order homogeneous linear equation over K is one of 
the groups in Sit's classification. 

Theorem 4.2. Let(DeC(t)(x)[8 x ] be an operator of order 1 or 2. Then, Gal A (CD/K) is defined over C(t). 

4.2 Improvements 

In computer algebra, one usually computes the partial-fraction decomposition of a rational function by com- 
puting the square-free factorization 13 of its denominator [10]. In the algorithms of §2, we stopped short of 
computing a complete square-free factorization of the denominators involved, because the polynomial d com- 
puted in each algorithm was sufficient to establish our results; namely, that the algorithms work, their outputs are 
defined over C(t), and that no factorizations of polynomials are necessary (beyond those which can be carried 
out using the Euclidean algorithm). The use of a single denominator d (and its powers) also allowed us to reduce 
the number of indices involved in our presentation. We believe that this has improved the clarity of our presen- 
tation. In practice, one should reduce the size of the system of linear equations considered in each algorithm by 
working with square-free factorizations of the denominators. 

Each algorithm in §2 constructs a first order inhomogeneous differential equation with undetermined co- 
efficients (these are the equations (H#), (Ijy), and (J,v)), and one then finds a minimal set of values for the 
coefficients such that the differential equation has a rational solution. One could also solve these differential 
equations abstractly, using any of the usual methods in computer algebra (e.g., Hermite reduction, Rothstein- 
Trager, Risch, . . .) [10], and then find a minimal set of values for the undetermined coefficients such that the 
solution is rational. Although any of these approaches would have to be abstractly equivalent to ours, they may 
well be more efficient in some situations. Therefore, the eventual implementation of our algorithms should first 
decide, perhaps heuristically, how to set up the systems of linear equations in order to minimize the amount of 
computation. 

4.3 Future directions 

The differential operators produced by the algorithms presented in §2 arise as solutions to a creative telescoping 
problem, namely: given rj £ F, where F is defined as in the beginning of §2, find L £ C(f)[3 t ] and / £ F 
such that L(r\) = 8 x f (see [4, §1] for a more precise and general definition); such an operator L is called a 
telescoper for T|. In [4] and [5], the authors propose algorithms which, in particular, solve the problem of finding 
a telescoper L £ C(f)[3 f ] for r\ £ C(x,t) arbitrary. Lemma 2.1 replaces a telescoping problem over F with a 
"twisted" telescoping problem over the smaller field K (see part (iii) of Lemma 2.1), for the restrictive class of 
algebraic functions r| £ C(x,t) such that rj" £ C(t)(x) for some n £ N, and also in the case that r| is transcendental 
over K, provided that ^ and ^ belong to C{t){x). It would be interesting to see whether one can use creative 
telescoping to compute PPV groups for higher-order systems, or over more general fields than we consider here. 
It would also be interesting to see whether our approach, which handles a very restrictive class of telescoping 
problems, can be generalized in order to "descend" some telescoping problems to simpler base fields, with 
Lemma 2.1 as a prototype. 

We mentioned in the introduction a result of [1 1], which says that the PPV groups we wish to compute are 
actually defined over C{t), while our algorithms find differential polynomial equations which are defined over 



13 lf F is a field of characteristic zero, and p £ F[x] is a polynomial, we say that p = ]T/ = i(P/) ! ls a square-free factorization of p if 
Pi £ F[x] for each (', and gcd t (p,,pj) = 1 for i 7^ j. The square-free factorization is unique up to multiplication by constant polynomials, 
and it can be computed using only the Euclidean algorithm and algebraic differentiation. 
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C(t). This is because we work with the solutions to 3. 1 found by Kovacic's algorithm, which are not (and cannot 
be) defined over C(t) in general. It would be desirable to have algorithms to compute the defining equations for 
the PPV group over the smaller field C(t), but we do not yet know whether it is possible to extend our methods 
to produce such an algorithm. 

It would be desirable to have algorithms in the case of several parametric derivations. We do not yet know 
whether it is possible to extend our methods to this more general setting in a straightforward way. In [12], 
the authors reduce the number of compatibility conditions that one has to check when working with several 
parametric derivations. It would be interesting to see if the results of [12] can be used to extend our methods to 
the setting of several parametric derivations, by working with one parametric derivation at a time. 
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